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NEW INVESTIGATION OF THE LA W OF ERRORS 
OF OBSERVATION. 



BY CHAS. H. KUMMEL, ASSISTANT TJ. S. LAKE SURVEY, DETROIT, MICH. 

1. Errors of observation being freed of constant and regular errors 
can be only subject to the laws of probability, and we may conceive an error 
A to be the result of two opposing influences ; one tending to make the ob- 
servation greater and the other smaller than the true result. If the increa- 
sing influence is the greatest a positive error is the result, and vice versa. 
Now conceive these influences to be made up from 28 = oo equal element 
errors ± i, some positive and some negative and let (8 -f v) times + i and 
(8 — v) times — i occur simultaneously to produce an error A then we have 

(1) (8 + v)i — (8 — v)i = 2vi = A. 

To produce an error A + dJ we must suppose one more + i and one less 
— i to preserve the constant total number 2i2, that is 

(2) (8 + v + l)i — (8 — v — l)i = 2(v + 1>" = A + dA. 

We have then 

(3) i — %dA, 

(4) vdA = A. 

We have here two events of equal probability, viz., + i and — i and 
their simple probability is = J. The terms of the development 

(«) s+r- 

will give then the probability of any number of + i and the remaining 
number of — i to occur simultaneously. The middle term, which is the 
largest, gives the probability of + 8i and — 8i occurring simultaneously 
or what is the same thing that of the error 0. If we denote this by <p and 
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so on the following terms on both sides by increasing integer subscripts, then 
<p v givis the probability of (Q + v)i and — (Q — v)i or by (1) that of the 
error A. We have then 

,„* -2« 2Q 2Q — 1 Q + l 

(6) f. =2 . T .- T - j- 

= probability of error = Qi — Qi, 
-2G 2Q 2Q — 1 Q + 2 

= probability of error dJ = (J2 + 1)* — (^ — 1)*> 

,~ -2a 2J2 2J2 — 1 i? + « + l 

(7) % = 2 . T .— g ^-^ 

= probability of error A = (£ + »)*' — (^ — «)*, 

(8) ?„ +1 -2 . T ..— g fl_„__i 

= probability of error A -\- dA = (Q -\- v -{■ l)i — (Q — v — l)i, 

<p _ = 2 _2fi . ^ = probability of error (Q—l)dA= (2Q—l)i—i, 

— 2S2 

(9) (p = 2 = probability of error i2c£ J = 2i2i — . i = <x>. 

These probabilities evidently form a continuous function. Denoting the 
general term <p v by f then <p H . 1 =<p + d<p, and dividing (8) by (7) we have 
probability of A + dA_ _ f v+1 __ <p + dxp _ i3 — v 
probability of A <p v <p Q + v -\- 1' 

df__ 2i?+l _ _ 2vdA + dA 
Ip Q + v + 1 X?dJ + «dJ +dzf 

««\ 2A + dA , ... 

( 10 ) = -M+TO3 by(4) ' 

and at the limit, since QdA = oo by (9), 

, , dg JiL = _ 2JdJ 

( Li) if ~ QdA QdA* 

QdA 2 is evidently a finite constant. Since QdA 2 = 4Qi 2 = 2 x 2Qi 2 , it is 
double the sum of the squares of the element errors ± i, and placing this 
quantity, viz., 

(12) 2Qi 2 = \QdA 2 = e 2 we have 

(13) *— ^ 
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Integrating we obtain log <p = log G — J 2 -5- 2s 2 . If J = then log <p 

= log (7, hence 

A 2 
(14) ? = ? e - 2^. 

We can determine <p from (6) for we have 

-2Q 2S 2£ — 1 S + 1 ,„ N 

-2a 1.2.3 2.G ,„ N 

= 2 - 12.2 2 .3 2 . ... £ 2 ( ^ = "^ 

By Wallis' theorem 

ff — 2 2 . 4 2 . 6 2 (2n — 2) 2 2w . _ . 

2 ~~ 1 . 3 2 . 5 2 (2» — l) 2 ' [n ' 

i su 2 . 4 . 6 . . . . 2ra , v 

{nn) = 1.3.5 ....(2n^ I) i (» " ~) 

(16) .-. p = (&r)r* (J2 = co) 

This remarkable result proves that the probability to commit no error at all 
is an absolute constant. 
We have by (12) 

therefore also 

(18) ^ = ^n) and 

/im dJ _A 2 

Let [j>]' denote the probability of an error to lie between a and b, then 

( 20) Ul = JJ7Wf ^ 

The sum of all the probabilities being certainty = 1, we must have 

/oi\ r ~i +co r +co ^ -4! ! 

and this is the case since by Laplace's integral 
(22) C"°dxe— X ' = i/n. 

J —00 
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2. Supposing now an ideal case in which all the possible errors, viz., 0, 
dA, 2dA, . . . . A, A -f- dA, . . . . (J2 — l)dA, QdA occur positive and nega- 
tive exactly in proportion to their probabilities as given by (19). Let m be 
the total number of errors then we have 

mdA _A 2 

m<p= WW) e Wi 

= number of times the error A should occur and 

A2 mdA ,<, _A 2 

= sum of the squares of the errors A. Following the usual Gaussian no- 
tation for sums of similarly formed quantities according to which, for inst., 

a i + a 2 + ■ • • • a m = [a 1 

a\ -\- a\ + al == [a 2 ], 

a 1 b 1 +a 2 b 2 -\- . . . a m b m — [a&], etc., 
we have the sum of the squares of all the errors in our ideal case 



[J 2 ] = mj - 7r77 ^e 2F 



_ 00 s 1 /(2n-) 

: f+ x 7 , A 2 

J JA.-e-^ 



me '' +eo 



~vW)l Je ^L„ + 7W>J d -i e * 



(-) - >=m- 



There can be no objection to use this formula for a finite number of er- 
rors provided we admit that in the long run the law of frequency of errors 
which we have assumed in this ideal case to be strictly fulfilled, will hold 
also in experience the more perfectly the greater the number of observations. 

As we shall see, the method which we are going to explain consists main- 
ly in this : to determine the unknown quantities in such a way that the sys- 
tem of resulting errors A 1} A 2 , . . . . A m resembles as much as possible the 
above ideal case, and this must bring us the nearer to truth the more regu- 
lar the law of error has been followed. 

The constant e found by (23) is called the mean error of A 1} J 2 , . . . A m 
and denoting it by e/± we have 
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and this quantity is to be distinguished from the mean error e of the obser- 
vations which slightly exceeds e/^ since it includes the uncertainty of the 
unknown quantities themselves. It is plain from all this that e the mean 
error of observation characterizes a class of observations, and this with re- 
gard to precision is plain from the following: If e' is the mean error of 
another system we have the probability of an error A' by (14) 

_A /2 
Y == f e 2f/2 - 
If an error A in the first system is to have the same probability as an 
error A' in the second then must 

A _A>___ 

£|/2 " e' 4/2' 
If then e' > e then also A' > J, that is, greater errors in the second system 
have the same probability as smller ones in the first. 
Gauss puts 

< 24 > h - TVS 

and calls h the measure of precision, which is the greater the lesser the class 

of observations. 

3. The most usual method to compare the precision of different systems 
is by means of the probable error. This is a quantity which stands in the 
middle of the series of errors arranged according to their magnitude. We 
have then the probability to commit an error less than the probable error 
equal to that to commit an error greater, or, denoting it by r, we have re- 
garding (20) and (21), 

/0 - r r dJ - A2 f" dJ - Aa 1 

(25) K^W) 6 ^^Kwwf "^4- 

Place 
(26) A-+-ei/2 = hA = d, 

then (25) becomes 

yrcJ V nJ 4' 

where we have put 
(28) r -5- ej/2 = p. 

If we solve (27) for p, which is most readily done by means of tables of 
the definite integral 

dJo e 


we find 



J 
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(29) p = 0.47694, 
hence by (28) 

(30) r = ep 4/2 = 0.6745s. 

This formula gives the relation betweed mean and probable error. We 
have then the following relation between mean error, probable error and 
precision 

(31) />=J_=£, 
and the probability of an error A is 

(32) f = -^e-^ = ^e- h2 ^ = <**<-£*' . 
v ' T e-|/(27r) j/n r\/it 

4. The probable error is highly significant in its relation to a determin- 
ed quantity. It is then an even wager that the error of that determination 
exceeds the probable error as that it is smaller. It is usually placed with a 
double sign after the qnantity. Thus the equation 

x 1 =a 1 ± r x 
signifies that the most probable value of the unknown quantity x x has been 
found a 1 with such an uncertainty that the actual error of a x is with the 
same probability greater as it is less than r v The probability of any val- 
ue to be = x x is therefore 

(33) fix,) = A-tr&b^)*. 
Similarly we have if x 2 = a 2 ± r 2 

(33') a(xa) = -£^a_ *-£(*»-«*»)*. 

If now 

(34) X = a x x x 4- a 2 x 2 

and we require the most probable value and probability of x we cannot as- 
sume the first to be a x a x + a 2 a 2> a * ^ eas * n0 * without proof. As this will 
be furnished as soon as we know the probability of X we shall determine 
this probability. 

"We evidently pass through all imaginable values for X if we combine in 

(34) any value of x x with any value of x 2 ; hence the probability of any 
special value of X is the compound probability 

p 2 dx x dx 2 < _jL( ai _ fll )s - £l(*,-a,)a 

after x x has passed from + 00 to — cxo, or 

(35) f (X) = ^^ f£e x g-&(*i-°x)*- 4(*.-«.) a . 
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By (34) we have 

(36) x 2 = X — aiXl 

,7V 

(37) dx 2 = 



*2 

dX 

since x 2 must be independent of x-^ after the x y integrations. Introducing 
these values into (35) we have 

(38) f{X) = - p2dX fdx 1 e-fj- (^-"i) 2 - 4?f (-T-«i*i-«.«».)". 

a i r \ T 1,' !Z —00 

The exponent may be put in the form — -4a;f + 25^ — (7, hence 
YK ' a 2 r x r 2 7tJ _J A A 

a 2 r x r 2 n A J _J A 

<W =^«- ( ^)-. > |i. [by (21).] 

But by comparison 

a =^ 2+ ^— <J~- = 7 i^[«?«M+(*-«.«»)'''! ] 

J g ^ngl ^(^f.^J^ Mfi -(X-« 2 a 2 )Vf "I 

With these values (39) becomes 

(40) f (X) = ^^2^-^fe^^---- ^- 

The most probable value of X is therefore 

(41) X = a 1 a 1 + a 2 a 2 , 
and its probable error 

(42) B = i/(«H + «!*•!, 
or the complete value 

(43) X = a x a t + a 2 a 2 ± l/(«i»1 + «2 r 2)- 
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If more generally 

(34') X = a 1 x 1 + « 2 a; 2 +. . . . a m x m = [ax], 

where a^ = a x ± r 1( x 2 = a 2 ± r 2 , . . . . * m = a m ± r m) 
then we have by composition 

< 40 ') ^>=^^-pw (X - [ " i,! > 

(43') Z=[ao]±T/([flV J ]). 

If 
(34") X^/^,^,....^) 

the above integration cannot be effected but an approximate solution can be 
given in that case which is the nearer perfect the smaller the probable er- 
rors r x , r 2) . . . r m . We have by Taylor's theorem, neglecting higher pow- 
ers of increments A#d l\x 2 , .... £±x m 
(34'") X = /K, a 2 ,...a m ) +f{a i )&x 1 + /'(a 9 )A* s + • • -/WA^ 

Within the range of A*i, A«2> A*V, for which this form is exact 

enough, X is of the form (34'). In the integration however these incre- 
ments have to pass from + co to — co. If the probable errors are small 
this will make no sensible difference since the integral 



f 



£^4.. 6 -4a 2 

A r\/ir 



approaches the more rapidly the smaller r. This circumstance admits to 
a certain extent the treating of X as a linear function of A«i, A#2> • • A«m 
and we have 

^4U j p^; l/([/'(«)V 2 » ^(°)' r *J 

(43") X =/( ttl , o a , . . . . a.) ± |/(L/V)V]). 

(lb 6e continued.) 



Note by S. W. Salmon. — In the note on Differential Calculus (p. 14), I 
wro te (y --A = 1. This needs to be proved. If the rate of mo- 

(y 1//\ • 
i— - — 7 J is greater 

than 1, and just after x = x', it is less than 1 ; therefore when x = x', 



(i 



— u 



-\ = 1. If B's rate is decreasing, it may be proved in a similar 



u — w J 
manner that [•" T ) = 1. 
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